356 


3H 3.2-9 


INVESTIGATION OF COSMIC KAI PROPAGATION 
IN INTERPLANETARY SPACE 

Chebakova E.A., Kolomeets E.V., Sevast'yanov V.N. 
Kazakh State University ,Timiryazeva St. 46, 
Alma-Ata 480121, USSR 


It was established experimentally that propagation of 
solar cosmic rays (scr) in interplanetary space up to 1 A.U. 
in most cases was of diffusion character ( Palmeira R.A.R., 
et al., 1971) . 

In present paper we describe solar cosmic ray events on 
the basis of equation given in the paper of Dolginov 
A.Z, ,Toptygin I,N.,1966)« 


an _ i 
W " 7 s 


fr»- ‘MB*** -^1} <» 


dr 


where n is density of particles in unit interval of kinetic 
energy £ ,r - heliocentric distance, E and P -total energy 
and momentum of particle correspondinly , V - speed of solar 
wind, 3B r - coefficient of radial diffusion, D* -coefficient 
of diffusion in space of energy ,(£=(2W,C 4 +£)/(fn,C >£)» where 
energy at rest. Equation (1) describes process of diffusion 
propagation of charged component of solar cosmic rays 
including convecton, adiabatic cooling and statistical 
acceleration. Numerical solution of equation (1) was based 
on the grid technique ( Krylov V.I.,et al.,1977 ). 


Boundary conditions were taken in the form: n(r,£ ,t)-0; 
CnV - 3B P || «Q(£,t) at r»r s , where Q is function of source, C 
-Compton-Getting factor, n(r. ,£,t)«0, r£ f r. ,r.l ; 

n(r ,£|, t}«=0, n{r ,£ 4 , t}~£”*® , £€[£*,£*] 

As a test problem we solved equation (1) by Monte Carlo 
technique using transit moments of the first and second 
orders of the equation (1) which is equation of Fokker- 
Planck type ( Jokipii J.R., Levy E.H., 1977). Figure 1 shows 
calculated time-intensity profiles at 11 MeV at r-1 A.U.. 

Solid line represents the results of calculation by grid 
technique, histogram-Monte Carlo technique for the case of 
instantaneous injection of solar cosmic rays at point n-0,01 
A.U. for y, =3 , 5 , r 4 -6 A.U. ,£<-0,2 MeV, £ 8 -50 MeV,#,,-^^, 
-const(r ),3e p (£-ll MeV)-l , 16x10** cm*/s, V-400 kra/s. 

Figure 2 shows calculated energy spectra If m obtained 
on the basis of maximum intensity at r-1 A.U. , when t-t m and 
histogram represents the results of Monte Carlo technique, 
solid line -calculations by grid technique, dotted line shows 
the form of energy spectrum at source. The solid line at' the 
bottom of figure 2 represents dependence of exponent J on £ 
,when Ifc m i3 approximated with the function of the form 
, dotted line is value of Jf® . Hardening of scr energy 
spectrum at 1 A.U. it^ dv^e to prevailing action of adiabatic 
cooling, when D £ «d ( Dolginov A.Z., Toptygin 

I.N.,1966), where AU-5tT kra/s - speed of alfven waves in solar 



wind plasma at r»l A.U.. 
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The dotted line in figure 1 is a result of solution of 
the equation (1) for all similar parameter and conditions 
except suggestion that 3Br<£)»const«l , 16x10** cn^/s. This 
condition is equvalent to % -const (t). Good agreement in 
behaviour of I(t) was obtained for t>t m when Jf »const(t ) 
f or aBp^const ( 6 ) . 

Usually when propagation of cosmic rays is investigated 
it is suggested that (r t £)mf(r) { f(£) , To study validity of 
this suggestion we described experimental solar cosmic ray 
time-intensity profile for 28,05,67 and 02,11,69 events. The 
data were taken from the paper of Lanzerotti L.J, ,1975, 
Description was carried out by minimization the difference 
between experimental value of I (t) and theoretical ones. 
Minimization was carried out for scr at given energies in 
suggestion of instantaneous injection of particles and. r t = 20 
A.U., when 3B P =3^r® at ranges 30,. ( r =» 1 A.U. ) 6 [3.75xl0 iS cn//s» 
3,75xl0 S1 # c^/s], b 6 [-2,2] with steps *«„= 0.1-**, Ab-0.5 , 

where <£ 0 -previous value of 8So, Figure 3 shows obtained by 
minimization dependences Aj»(£), b(£) for 28,05.67 event and 
dependences A P (R), b(R) for 02,11.69 event, where A p -radial 
mean free path, R- rigidity of particle ( A P *33B r /l/ , IT 
-particle speed ). 

Figure 3a shows values of A P (£) and 6(6) obtained by 
minimization of Ot*(circles) and by least square method 
(rectangles) for Jf 0 “3 for protons at energy £ : 1.2-2, 4 

MeV,2.5-4.3 MeV, 4. 4-5.0 MeV; 5. 0-9. 4 MeV; 9.4-17.4 MeV; 

16.5- 19.7 MeV. Figure 3b shows ?| P (R) and b(R) for 02.11,69 

obtained by mininraization using least square technique for 
electrons at energies: >0.35 MeV; >0.6 Mev; >1,1 MeV; for 
protons at energies: 1,1-2, 5 MeV; 2, 5-4, 3 MeV; 5. 0-8, 8 MeV; 
5. 9-8. 8 MeV, 8.8-16,7 MeV, 17.0-19,7 MeV; for alfa-particles 
at energies: 3.8-6. 2 MeV; 6. 2-8. 5 MeV, 8.5-17.5 MeV, 

17.5- 24,5 MeV, 24.5-42.5 MeV, 42.5-83.5 MeV. Rhombs, circles 
and rectangles in the upper part of figure 3b are values of 

A P for electrons, protons and alfa-particles correspondingly. 

Analysis of figure 3 reveals that h r ~R* , 4*0,6 , 

4*0,27 for events 28.05.67 and 02.11.69 correspondingly and 
98^R>0. Solid smooth curves in figure 4 show the results of 
description of time-intensity pr.ofile of scr protons for 
28.05.67 for Jfo-2.3, « P = 3 e o r*£ , 4=0.8, b=0. 75 ^6 , 

[6] «* MeV, Coordinate axis for the given curve is shown 

with the help of arrow. Figures near curves are energy 
ranges of registration of solar cosmic rays. Solid line in 
figure 5 is description of energy spectrum of proton 
intensity If (8)for the given event. 

Prenounced softening of It m (8) at r*lA.0, 

(If ) with respect to scr spectrum at source where 

Kr*0;r»r $ )~d’ i * >should be noted. 
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FIGURE 5 
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Analysis of abovementioned results shows that 
suggestion of instantaneous injection allows to describe 
propagation of scr when function SBp C r, R) can not be 
factorized with respect to r and R and of Xp^P'R ,then 
d$/SR>0 and at r«l A.U, <4«0.3-0.6 that is in good agreement 
with the data on frequency power spectrum of interplanetary 
magnetic field fluctuation. 
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